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Abstract. - We propose a cooling scheme based on depolarisation of a polarised cloud of 
trapped atoms. Similar to adiabatic demagnetisation, we suggest to use the coupling between 
the internal spin reservoir of the cloud and the external kinetic reservoir via dipolar relaxation 
to reduce the temperature of the cloud. By optical pumping one can cool the spin reservoir 
\£) , and force the cooling process. In case of a trapped gas of dipolar chromium atoms, we show 

. that this cooling technique can be performed continuously and used to approach the critical 

f*"*) ' phase space density for BEC. 
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o 
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Introduction. - Adiabatic demagnetisation [1,2] is a well established and very efficient 
Q-i, cooling scheme which enables researchers in solid state physics to cool their samples by several 
orders of magnitude in a single cooling step [3,4]. However, depolarisation processes have 
not yet led to a cooling concept in atomic physics. Instead, evaporative cooling which can be 
observed in many fields of physics is typically applied to obtain temperatures in the nK regime. 
This cooling mechanism was proposed and demonstrated for magnetically trapped atoms by 
Hess [5]. Meanwhile, it has been studied intensively [6,7] and could also successfully be applied 
to atoms [8,9] and molecules trapped in optical dipole traps [10]. Up to now, this scheme 
is essential to obtain degenerate atomic or molecular quantum gases and allows nowadays to 
generate gases with temperatures below T = 500 pK [11]. By means of a controllable finite 
trap depth Uq, high energetic particles carrying more than the mean energy of a trapped 
particle are allowed to escape from the trap. Rethermalisation of the remaining particles via 
elastic collisions reduces the temperature of the trapped atomic cloud and at the same time 
produces particles which have sufficient energy to leave the trap again. In typical experiments, 
this technique allows one to increase the phase space density p — no(2Trh 2 /(niksT)) 3 / 2 by 
several orders of magnitude, where no and m denote the peak density and the atomic mass, 
respectively. The ratio r\ ev — Uo/iksT) between the trap depth and the thermal energy of 
the cloud is commonly referred as cutoff parameter. The higher this ratio is chosen, the 
more energy can be carried away by a single atom and the less atoms are lost to achieve 
the final temperature. However, in this case the particles need more time to rethermalise, 
so that trap losses become more significant and finally limit the cooling process. Thus, the 
efficiency x 01 the cooling process is defined by the gain in phase space density per atom loss 
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(— d(ln p) / d(\n N)) which can be optimised during the evaporation using the cutoff parameter. 
Typical experimental values of x range up to 4. Beside the intrinsic high loss of atoms of this 
cooling method, in optical traps the trapping volume has to be enlarged in order to reduce 
depth of the trapping potential, thus the forced evaporation cooling typically does not reach 
the runaway regime where the evaporation accelerates itself. 

To reduce the temperature of an atomic sample more efficiently, we suggest to transfer 
kinetic energy to an internal degree of freedom (spin) while the cloud depolarises via inelastic 
dipolar relaxation collisions. Subsequent optical pumping connects the spin reservoir to the 
light field and allows to remove the energy from the trapped atomic cloud. This cooling 
scheme does not rely on removing atoms from the sample or changing the trapping potential. 
Therefore, it is expected to be much more efficient than conventional evaporative cooling. 

In the following, we consider N atoms in an homogenous magnetic offset field B trapped in 
a power-law potential, which is independently of the internal state characterised by U(x, y, z) = 
c x x ni + Cyy 712 + c z z n3 with a = nj 1 and may be realised by a far off resonant optical 
dipole trap. Having in mind a specific element - 52 Cr, which has been Bose-Einstein con- 
densed recently [12] - we focus in this letter on atoms in a stable state without both hypcrfinc 
structure (/ = 0) and electron orbital momentum (L = 0), though these are no restrictions 
for the cooling scheme. A finite electron spin 5 leads to 25 + 1 magnetic substates \mg) 
which are energetically separated by the Zeeman energy AE Z = 2[i B B (see fig. P, where \i B 
is the Bohr magneton. The dipole moment can cause inelastic dipolar relaxation collisions 
in which the total spin quantum number of both atoms is not conserved. Energy conserva- 
tion requires, that for each spin flip event to a neighbouring lower (Arris — ~ 1) or higher 
(Arris — 1) energetic state the energy AEz is transferred to or detracted from the kinetic 
energy of the colliding atoms. Moreover, we assume for the simulation, that thermalisation 
between the external degrees of freedom (kinetic reservoir) via elastic collisions occurs much 
faster than the energy transfer to the spin reservoir via dipolar relaxation. Such we can im- 
ply thermal equilibrium of the kinetic reservoir and the total energy of the kinetic reservoir 
is given by E = (3/2 + a)Nk B T. Starting from a non-equilibrium distribution of atoms 
across the states \ms), the sample gradually relaxes via dipolar relaxation to equilibrium 
occupation of the states which is given by the Boltzmann distribution. For a net relaxation 
rate N r = 53-t=i ^S,* which contains spin flip events to neighbouring states, the cooling rate 
can then be estimated using the time derivative of the total energy of the kinetic reservoir 
(E = (3/2 + a)k B (TN + TN)): 

AE z N r 

T= (1) 

(3/2 + a)Nk B ' 1 ] 

where we have neglected atom loss. Depending on the starting condition relaxation will cool 
or heat the sample. The latter process limited us to obtain Bose-Einstein condensation in a 
cloud of chromium atoms in a magnetic trap [14]. 

If the atoms with temperature T are polarised in the energetically lowest state, transitions 
to the neighbouring higher energetic substate caused by dipolar relaxation collisions cool the 
sample. As the sample approaches equilibrium, the relaxation rate and therefore the cooling 
rate tend to zero. The final reachable temperature T eq for a certain magnetic field is depicted 
in fig. and can be calculated from: 

6 I „, 1 Z,_T1 _ AT I I I „, > ;„_T I AZ7„^i=0 e 



E = N[-+a\k B T =N\[-+a\ k B T eq + AE Z ^=» A — | . (2) 
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Fig. 1 - Dipolar relaxation collisions into energetic higher state. Shown are the lowest two energy 
levels of a ground and excited state manifold of a J — > J' = J-transition, respectively. The levels 
are separated by the Zeeman energy AEz- Also indicated is the distribution function f(E) of the 
relative kinetic energy of a trapped polarised atom cloud in thermal equilibrium. Inelastic single 
(ssf) and double (dsf) spin flip transitions are only possible for atoms in the high energy tail of the 
distribution. Using <r~-polarized light the sample can be polarised in the energetically lowest and 
dark state (|ms = —S)). 

Fig. 2 - Equilibrium temperature of a previously polarised cloud depending on the magnetic field 
and the spin quantum number S. The initial (a) and equilibrium (b) situation of sample with S=3 is 
indicated in the inset. 



Unlike in solids used for adiabatic demagnetisation where the phonon heat capacitance in 
a cryogenic surrounding (To ~ IK) is negligible compared to the heat capacitance of the 
spin reservoir, in a trapped atomic gas (To ~ 1 — 10 3 /i_R") the heat capacitances of the spin 
and kinetic reservoir are of comparable magnitude. Therefore, the achievable temperature 
reduction is much smaller and an optimum (T eq /To = (3/2 + a)/(5/2 + a)) is theoretically 
obtained for S — ► oo in the limes of B — > 0. However, in atomic physics the spin reservoir 
can be cooled very easily by optical pumping. In this way, this depolarisation process can be 
repeated several times or even driven in a continuous way, like we will discuss in the following. 
Experimentally, polarisation of the sample in the energetically lowest state \ms) = —S can 
e.g. be accomplished using er~-polarised light on a J = S - > J' — J transition. In this case, 
| rris) = — S is a dark state and its population is not affected by the pumping light (see fig-QJ. 

In a cycle consisting of a dipolar relaxation collision and an optical pumping transition, 
cooling can be provided if the Zeeman energy AEz exceeds the energy (E po i ~ E rec = 
k 2 h 2 /(2m)) needed to polarise the cloud. In principle this scheme allows to generate samples 
with temperatures below T po ; = E po i/kB, since in a thermal distribution high energetic atoms 
E > AEz > Ep i will undergo dipolar relaxation collisions and contribute to the cooling of 
the cloud. In this respect, dipolar relaxation can be considered as an evaporation out of the 
energetically lowest state, where the cutoff parameter is given by r/s — AEz / (ksT). Note, the 
energy reduction of the sample works without thermalising collisions. The required ingredients 
are a high inelastic dipolar relaxation rate and a optical pumping transition back to the initial 
state, which is a dark state for the pumping light. A scheme based on elastic collisions (spin 
changing collisions) in combination with a quadratic zeeman shift of the magnetic substates 
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was proposed by G. Ferrari [13]. In his case, either linear or circular polarised pumping light 
provides the necessary pumping mechanism including a dark state. 

Cooling model. In the following, we develop a model to describe the continuously 
driven cooling process and study the practicability of the scheme for a sample of chromium 
atoms. Dipolar relaxation rates which will be used in this section have been previously studied 
experimentally and theoretically in our group [14]. There we found that the process is well 
described by dipole-dipole scattering in the first Born approximation where either no (elastic 
collision) , one or both colliding atoms undergo a transition to a neighbouring substate (Amg = 
0, ±1). If AMg — Anig i+Ams.2 = —2, 2 accounts for the total change in the spin quantum 
number of both atoms, AM$ AEz is the energy which is released (sign(AMs) < 0) or required 
(sign(AMs) > 0) during such a collision. The inelastic dipolar relaxation rate for an atom 
cloud polarised in a extremal magnetic substate into the neighbouring state is given by N r = 
Ndip — —fiN 2 /V , where (3 = {(u\ + 2o~2)v re i) therm is the thermally averaged rate constant 
containing single and double spin-flip transitions and V is the mean trapping volume. Hereby 
we introduced the collision cross sections for single spin-flip (<ri = £5 3 (1 + h(kf t i/ki))kf t i/ki) 
and double spin- flip (o- 2 = ^S 2 {l + h{k f ^/ki))k f ^/ki) events with £ = (/zo(2^ i3 ) 2 TO) 2 /(307rft 4 ) 
and the relative velocity v re i, respectively. h(x) includes the symmetries of the particles and is 
defined in [14]. Finally, the factor kf t AM s /ki accounts for the different density of final states 
in the inelastic process and therefore ensures that the energy conservation is fulfilled: 



!.„ A „ f /i in AM s AE Z :f i rnAMs AEz 

«/,AM S = I y 1 ! pfcf — " L> — plf — ' (3) 

ki y otherwise. 

Thus, the energy transferred to spin reservoir via dipolar relaxation collisions is given by 
Edip — s\gTi(AMs)AEzNdi P - While AEz linearly increases with the magnetic field, Ndip 
decreases for atoms polarised in the energetic lowest substate (\ms) = —S). Therefore, the 
amount of transferred energy can be optimised by the cut-off parameter r]B- In the case of 
trapped chromium atoms (S=3), it is maximum for r)B,opt ~ 1-31. If the temperature of the 
cloud exceeds the recoil temperature (T 3> T rec ), other heating mechanisms can be neglected 
and we obtain a rough estimation ( 1 ) for the expect cooling rate: 

T(vb) « -^—J^tS 2 {(1 + Vb)S + (2 + 4r, B )e-^ } VB e~^ =T 3 / 2 . (4) 
3/ 2 + a V 7rm L ' V 

Since the rate intrinsically depends on the dipolar relaxation rate, the cooling process works 
especially well for dense samples of atoms or molecules which possess high spin quantum 
numbers. The sum within the curly brackets contains both spin flip transitions whereas 
the contribution of the double spin flip transition is suppressed by a factor of e~ VB . More- 
over, we find that the time evolution of the temperature is closely connected via the mean 
volume (V oc T a ) to the form of the trapping potential since a reduction in temperature 
simultaneously leads to an increase in density. While a three dimensional harmonic potential 
(V = ( V 'AwksT I (ZJ^/m)) 3 , where ZU is the mean trapping frequency) results in a linear change 
of the temperature with the volume, a potential with a > 3/2 - like it can be realised in a 
dimple trap - enhances the cooling process in time and yield a runaway behaviour. 

In order to obtain a more realistic model, we include atom loss and the optical pumping 
process. Hcrby we will restrict our following considerations to a three dimensional harmonic 



( 1 )To calculate r)g^ opt and eq. J1J we approximated (l+h(kf t 2/ki))kf t &M s /ki = l/2-e(«„rv' 2 fB AM sH 
for atoms in the lowest state. 
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trapping potential. Atom loss processes caused by background gas collisions and three-body 
recombination are characterised by the rate constants l/rhg and L 3b , respectively. We do not 
expect two-body losses, if a sufficient deep potential is provided (Uq ^> AEz and N ev <C Ndi P , 
where N ev is the evaporation rate out of the trap). The evolution of the total atom number 
in the trap is then described by: 

N = N bg + N 3h = - — N-^§N 3 . (5) 

Ug V 

Depending on the loss process the kinetic energy of the cloud is changed on an average by: 

Ei oss = ik B TN bg + 2k B TN 3b . (6) 

Since atoms lost in a three-body collision carry less than the mean energy of a trapped gas, the 
remaining cloud heats up. To account for polarisation effects, we consider the two energetic 
lowest states (\ms — —S), \ms = —S + 1)) and assume a negligible population in the other 
states if the cloud is optically pumped. In the following, indices 1 and 2 will indicate these 
states, respectively. The net exchange of atoms due to dipolar relaxation is given by: 

N r = N d ip,1^2 ~ Ndip,2^1 (7) 

where Ndip,i—»2 & n d Ndip,2—fi are dipolar relaxation rates containing collisions between the 
atoms in the same state and atoms in state 1 and 2 which result in increase or decrease of 
atoms in state 1, respectively [15]. Including the pumping process, the population in state 1 
and 2 evolve according to the following rate equations: 

N 1= N r + Nj± + ((1 - k)N 2 - pN^Tsc, (8) 

N 2 = -N r + - ((1 - K )N 2 - pN^Tsc (9) 

Hereby, we introduced the scattering rate r sc and the Clebsch-Gordan coefficient for a 
er~ -transition from state 2. The term pNiT sc takes into account that in the experiment the 
light is not perfectly polarised. Absorption and emission of light are each accompanied by a 
momentum transfer. Since the light scattering rate is on the order of the dipolar relaxation 
rate which is according to our assumptions much smaller than the elastic collision rate, the 
sample thermalise between two scattering events and the net energy transfer of an absorption 
spontaneous emission cycle reads: 

E P oi = {pNx + N 2 )E rec T sc . (10) 

These equations neither consider effects of quantum degeneracy nor reabsorption of scattered 
light. In particular the latter will limit the cooling process of an extremely dense sample. 
However, effects of reabsorption can be reduced in a lower dimensional trap or in the so called 
festina lente regime [16], where the trap frequencies exceed the linewidth of the pumping 
transition. Thus, the change in the kinetic and potential energy of the atomic cloud is given 
by: 

E = E dip + Ei oss + E po i. (11) 

We solved the equations iJSJ and 111( 1 numerically and optimised ijb after each time step 
to achieve maximum efficiency x- The calculated time evolution of the required magnetic 
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Fig. 3 - Time evolution of temperature (solid line) and magnetic field (dashed line) during the 
simulated cooling process. Additionally indicated are the solution of eq. 2] (dotted line) and the 
temperature (dashed-dotted line) corresponding to the energy needed to pump atoms from state 2 
into state f by <r~-polarised light. The same data are shown in the inset on a different time scale. 



field and the temperature of a sample of 5 • 10 6 chromium atoms with an initial temperature 
of 200 fiK in a harmonic trapping potential (ZU = 2tt ■ 500 Hz) are depicted in fig. Hereby, 
we assumed p = 10 -3 and altered the scattering rate T sc between 30 — 2000 s _1 to fix the 
population in state 2 to 2% of state 1 during the simulation. Since additional heating caused 
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Fig. 4 - Atom number and density (inset) dynamics during the cooling process. The solid line 
correspond to the curves shown in fig-OKp = 1CT 3 , N 2 /N 1 = 0.02). The dotted line and the dashed 
curves represent simulations for p — 10~ 2 which accounts for the polarisation of the pumping beam 
and N2/N1 = 0.005, respectively, while keeping others parameters fixed. 

Fig. 5 - Phase space density gain per atom loss. The displayed curves correspond to the curves 
illustrated in fig.01 The efficiency \ given by the slope of the curve, reaches a maximum value of 570 
for the solid curve. Also indicated is the area (gray) where typical evaporation cooling trajectories 
would be expected. 
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by atom loss and optical pumping are negligible during the first 7 s, both the magnetic field 
(dashed line) and the temperature (solid line) reduce linearly in time (inset fig.[3J). This process 
is well approximated by eq. which is indicated by a dotted line. During this time, atom 
loss is mainly given by background gas collisions (rtg = 200 s) (see fig. El solid line). The linear 
decrease in temperature leads to diverging increase in both peak density (oc T -3 / 2 ) (see inset 
fig. 131 (solid line)) and phase space density (oc T~ 3 ), until three-body recombination finally 
limits the achievable density of the atomic cloud. Preliminary experimental measurements on 
a trapped cloud of Cr atoms suggest a rate constant L^t, ~ 1 0~ 41 m 6 , which prevents us from 
obtaining densities much higher than 10 15 atoms m~ 3 in the simulation. A clear signature for 
these collisions is the more pronounced atom loss which can be observed after 7 s. 

The density still rises for a little while, until heating processes have the same magnitude 
than the cooling processes. The cooling gets inefficient and a final almost constant temperature 
close to the recoil limit is reached at a magnetic field of about 10 mG. In our specific case the 
final temperature is below the temperature which corresponds to the mean pumping energy 
(Epoi = k3 Erec) and which is indicated by a straight line in fig. 

In fig. [o] we illustrate the cooling process in a double logarithmic phase space density - 
atom number plot (solid line). The slope represents the efficiency of the cooling process as 
defined before. Until three-body collisions and heating caused by optical pumping limit the 
process, 5 orders of magnitude in the phase space density could be gained while only 5 percent 
of the atoms are lost. Here we obtain values of \ starting from 90 to about 570. Then the 
efficiency drops and finally the phase space density decreases again. 

Fig. ^ and fig. [S] contain additional curves which result from simulations in which we 
changed either the degree of polarisation of the pumping beam p = 10~ 2 (dotted line) or the 
ratio N2/N1 = 0.005 (dashed line). In the first case the pumping process is worsened and 
additional photons have to be scattered to obtain the ratio N2/N1 = 0.02. The extra heat 
reduces the cooling rate, so that density increases slower and atom loss is less significant. 
The final reachable phase space density is reduced by 1.5 orders of magnitude. Due to the 
improved polarisation of the sample, the cooling rate is increased in the second case and the 
density rises faster in the beginning of the process . However, since more photons are required 
to maintain the degree of polarisation of the sample the final temperature (~ 1.7 fiK) is higher 
and the corresponding density is lower. 

Conclusion. We studied a cooling scheme which is based on depolarisation of a po- 
larised cloud via dipolar relaxation collisions. Our results imply that for a cloud of optically 
trapped chromium atoms with experimentally reachable starting conditions and parameters, 
the cooling scheme works extremely efficient and may allow to generate a degenerate sample. 
Even if this goal can not be achieved, excellent starting condition for evaporative cooling are 
obtained. The scheme may also be applied to other atomic (e.g. He, Cs) or molecular species 
which exhibit large inelastic dipolar relaxation collision rates between magnetic substates and 
which can be polarised optically. Especially for molecular samples it might be a way to in- 
crease the density be serval orders of magnitude. The scheme allows to remove energy from 
a trapped sample and to increase its density without the need of thermalising collisions. The 
scheme may also be interesting for sympathetic cooling of mixtures of species (e.g. Cr-Rb [17]). 
Interspecies collisions may transfer energy and spin, so that optical pumping has not to be 
applied to the species which undergoes dipolar relaxation. 

* * * 
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